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Banderillas: Effects of deforestation on dragonflies 
(Insecta, Odonata) in the Pacific lowland of Costa Rica 

Banderillas: Efectos de la deforestacion sobre libelulas 
(Insecta, Odonata) en la tierra baja Pacifica en Costa Rica 

Florian P. Hofhansl & Stefan Schneeweihs 


Abstract: Deforestation in the past is responsible for the loss of large forest areas in all tropical regions and this process is ongoing. 
The reclamation of cultivated land reinforces the disappearance of natural habitats and reproductive sites for many species typi- 
cal for pristine forests. Lotic waters may remain when the trees are gone, but their characteristics will have changed dramatical' 
ly and species tied to these water bodies face different challenges maintaining self-sustaining populations. This study contributes 
to the question of whether changes in land use affect species assemblages and their ability for recruitment, and whether streams 
embedded in cultivated landscapes provide viable reproductive habitats for dragonflies. 

We surveyed dragonflies at four forest sites, two sites at forest margins and four sites at streams in an agricultural landscape. At 
each site, adult dragonflies were recorded twice along a 50 m transect in January and February 2007. We found significant diffe¬ 
rences in species richness, species assemblages and proportion and number of widespread species between forest and cultivated 
area. Only 3 of 11 species found in forests were also observed at streams in agricultural areas. Our results indicate that rivers in 
agricultural environment do not represent suitable reproductive habitats for the majority of ‘forest species’. Therefore, deforesta' 
tion poses an enormous threat for these species. 

Key words: biodiversity, conservation, landmse, species composition, species richness. 

Resumen: En Costa Rica y otros pafses tropicales, la deforestacion es responsable de la perdida de grandes areas de bosque que son 
habitats naturales de varias especies. Con la explotacion de madera y el cultivo de areas antes forestadas van desapareciendo los 
lugares donde viven y se reproducen estas especies. En el caso de las libelulas los habitats de reproduccion se encuentran cerca y 
dentro del agua, donde tambien viven sus larvas. Los adultos se encuentran generalmente en las orillas para aparearse y finalmen' 
te las hembras ponen sus huevos cerca del agua o directamente dentro. A pesar de la deforestacion aun permanecen los arroyos 
en las zonas cultivadas. Sin embargo sus caracteristicas cambian y las libelulas enfrentan situaciones totalmente diferentes. 

Se investigaron las diferencias entre las comunidades de libelulas de arroyos del bosque y zonas cultivadas, y si las aguas en zonas 
cultivadas sirven como lugares de reproduccion para libelulas. 

Se observaron libelulas adultas en cuatro puntos de un arroyo en el bosque, en dos arroyos situados al linde del bosque y en cua' 
tro arroyos en la zona cultivada. Se visitaron cada sitio dos veces a lo largo de un transecto de 50 metros, durante enero y febre' 
ro de 2007. La composicion de libelulas variaba significativamente entre el bosque y la zona cultivada, en la cual habia mas db 
versidad de especies. En la zona cultivada habia tambien un porcentaje mas alto de especies de distribucion amplia y un mayor 
numero de individuos de esas especies. Solo 3 de las 11 especies encontradas en el bosque se encontraron tambien en los arroyos 
de la zona cultivada. Los resultados senalan que las especies sensibles del bosque no pueden utilizar las aguas en zonas cultivadas 
para su reproduccion y que por tanto la deforestacion representa una gran amenaza para estas especies. 

Palabras clave: biodiversidad, conservacion, uso del suelo, composicion especifica, riqueza de especies. 


Introduction 

Banderillas, the local name for dragonflies (includ' 
ing damselflies) in southern Costa Rica, are among the 
most common and attractive groups of insects. Their 
aerial acrobatics, beautiful colours and relatively large 
size capture everyone’s attention. Despite a widespread 
ignorance about these creatures, they might easily be 
used as indicators to quantify effects of changing envi¬ 


ronmental parameters in freshwater habitats (e.g. oxy' 
gen concentration, nutrient accommodation and pollu¬ 
tion), since dragonflies in all stages of their life cycle 
are closely bond to aquatic habitats (ESQUIVEL 2006). 

Odonata inhabit a variety of habitat types: larvae 
are capable to develop along the whole course of fluvial 
systems, from headwaters to estuaries, and in lentic wa¬ 
ter bodies, both perennial and temporary (CORBET 
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1999). However, among dragonflies habitat require¬ 
ments can vary significantly. Some species show prefer¬ 
ences for a certain habitat type and can only persist 
within narrow niches of biotic and physical factors. In 
contrast, other species are eurytopic and therefore rela¬ 
tively tolerant to a broad range of environmental vari¬ 
ables and thus inhabit different types of freshwater habi¬ 
tats. Usually, the physiological and morphological con¬ 
straints of the larvae determine the presence of a species 
at a given location. Nonetheless, it is the adults that 
oviposit and their habitat selection acts as a filter on 
species composition (McCauley 2006, Hofmann & 
MASON 2005). However, reproductive and foraging 
habitats of dragonflies are affected by human activities 
and changes in land use. In Costa Rica, like in most 
tropical countries, deforestation converts forested areas 
into open and fragmented landscapes (SANCHEZ- 
AZOFEIFA 2001). When forest is converted into agricul¬ 
tural landscape, the characteristics of terrestrial and 
aquatic habitats (e.g. flow modification, nutrient and 
pesticides influx, sediment deposition, invasion of alien 
species in riparian vegetation and isolation due to land¬ 
scape fragmentation) change (DUDGEON et al. 2005). 
Such severe alterations in environmental parameters 
are likely to have a remarkable impact on dragonflies 
and their occurrence in modified landscapes. 

Deforestation and fragmentation obviously have a 
major impact on ecological systems. The composition 
and configuration of countryside habitats strongly in¬ 
fluence the diversity and composition of native plant 
and animal communities (SWEENEY et al. 2004, DAILY 
et al. 2003). The consequences of forest modification 
and habitat conversion on species richness have been 
the subject of intensive research at a multi-taxon level 
(Lawton et al. 1998, Schulze et al. 2004). In forests, 
slight disturbance may increase species richness, but 
can also cause a shift in species assemblage and extinc¬ 
tion of sensitive species (BROWN 1997). However, in 
strongly modified landscapes, species diversity usually 
declines significantly and sensitive forest species are re¬ 
placed by a few generalists and geographically wide¬ 
spread species (SlGEL et al. 2006). For tropical (CLAUS- 
NITZER 2003, SAMWAYS 2003) and subtropical (Stew¬ 
art & SAMWAYS 1998) regions, studies on dragonflies 
reported similar results. In Costa Rica, studies have fo¬ 
cused on mammals (DAILY et al. 2003), birds (DAILY et 
al. 2001, Matlock et al. 2002), butterflies (Horner- 
Devine et al. 2003), moths (RICKETTS et al. 2001) and 
other arthropod groups (GOEHRING et al. 2002), but 
knowledge of dragonflies is scarce. In this preliminary 
study, we present initial data on the response of dragon¬ 
fly communities to forest modification and conversion 
in Costa Rica. 


Particularly the following life-history traits of 
Odonata may be responsible for their sensitivity against 
habitat disturbance: 

(1) Species often depend on certain types of macro¬ 
phyte vegetation or even single plant species 

(Schindler et al. 2003, Rantala et al. 2004). 

(2) Thermoregulatory responses to open or closed 
and shaded habitats may differ between the two subor¬ 
ders. Smaller Zygoptera have lesser thermoregulatory ef¬ 
ficiency than Anisoptera which can adjust their heat 
balance by wing motion or body posture in relation to 
the sun (De Marco et al. 2005). 

This study addresses the question of how much a 
changing environment modifies dragonfly assemblages 
in the Pacific lowlands of Costa Rica. Our assumptions 
are that sensitive semi aquatic insects like Odonata di¬ 
vide into species that cannot live or reproduce in dis¬ 
turbed habitats, and species overruling previously dom¬ 
inant species in abundance in these new, artificially 
shaped niches. 

Material and Methods 

Study Area 

Dragonfly surveys were conducted in the vicinity of 
the Tropical Research Station La Gamba (8°42’46”N, 
83°12’09”W) in the south of Costa Rica, in the ‘Parque 
Nacional Piedras Blancas’, a protected area in the Golfo 
Dulce region comprising one of the last remnants of low¬ 
land rainforest on Costa Rica’s Pacific coast, and adjacent 
cultivated areas. Average annual precipitation in the re¬ 
gion is 5000-6000 mm with no month receiving less than 
180 mm. The rainy season lasts from September until 
November, the driest months are between January and 
March. Average monthly temperature ranges between 
22°C and 28°C (Weber et al. 2001). The Golfo Dulce re¬ 
gion with its high richness of vascular plants (e.g. 700 
tree species) is a biodiversity ‘hot spot’ and a high degree 
of endemism corresponds with the region’s status as a 
refuge during glacial periods (WEBER et al. 2001). Beyond 
the borders of the national park, human impact on the 
environment becomes obvious as the landscape changes 
dramatically from forest to a matrix of agricultural areas 
consisting of cultivated fields, plantations and pastures. 
The main cultivated plants are rice (Oryza sp.) and oil 
palm (Elaeis guineensis ). In these cultivated areas, trees 
have almost completely vanished and often only gallery 
forests of several metres width remained, accompanying 
bigger streams as ‘forest like’ structures. 

Study sites 

To investigate differences in species assemblages be¬ 
tween pristine and modified habitats, we surveyed drag- 
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Fig. 1 : Position of 
rivers and survey sites. 
Dark areas are 
covered with forest, 
pale areas are open 
landscapes. 


onflies along streams at four pristine forest (fl-4), two 
forest margin sites (fml-2) and four sites within the cul¬ 
tivated area (cal-4) (Fig. 1). 

Pristine forest sites were situated along first to sec¬ 
ond order streams. Three sites were located along the 
river ‘Quebrada Negra’ (Fig. 9) and one at a small 
stream in a narrow and shaded ravine. 

The two forest margin sites were situated at streams 
close to the garden of the field station. Fm2 at the Que¬ 
brada Negra was surrounded by the garden of the field 
station and an old agroforest with Theobroma cacao . Tall 
trees (e.g. Tectona grandis ) provided a high degree of 
shading. In contrast, fml was located at a small stream 
supplying an artificial reservoir. This stream runs 
through a meadow that is occasionally mown, without 
any trees in the riparian vegetation. 

Cultivated area sites were located along four differ¬ 
ent streams: Rio La Gamba, Rio Bolsa, Rio Chorro & 
Rio Oro in the hinterland of La Gamba village, running 
through flat, human altered landscape (Fig. 8). Riparian 
vegetation consisted of gallery forest strips with tall 
shading trees (e.g. Luehea seemanii & Inga sp.) and vari¬ 
ous macrophytes covering the banks. 


Survey Methods 

Between January 16 th and February 17 th 2007, each 
site was surveyed two times by pacing the riverbed for 
one hour within 50 metres of the shoreline. Dragonflies, 
either flying by or perching on vegetation at both sides 
of the stream, were caught with a sweep net (mesh- 
width 3 mm, diameter 45 cm) fixed on a stick (1.5 m). 

Dragonflies were identified according to FORSTER 
(1999) and afterwards released to their former territo¬ 
ries. Identification was always possible to genus level 
but sometimes not to species level. In such cases, we as¬ 
signed individuals to ‘morphospecies’ (KRELL 2004). 
Carlos Esquivel (Universidad Nacional Costa Rica, 
Heredia) kindly provided help in species identification. 

Statistical analysis 

Species richness 

To assess species inventories of habitats we conduct¬ 
ed observation units. These observation intervals are 
temporary and therefore only represent a snapshot of 
species occurring at a given location. Results are strong¬ 
ly influenced by the mobility of species, external factors 
(e.g. weather conditions) and stochastic events (BROSE 
& Martinez 2004, Rouquette & Thompson 2007) and 
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Fig. 2 : Box plots of 
recorded species in 
forest (f: n = 4), forest 
margin (fm: n = 2), 
and cultivated area 

sites (ca: n = 4) co 10 
habitats. Dots show q 
median species q_ 
number per habitat 8 

type, boxes show 25- 
75% of recorded o 
species per habitat ^ 
type and whiskers 6 

show minima and 
maxima of species 
recorded at sites 4 

within one habitat 
type. 

it is very likely that they underestimate species richness 
at the surveyed sites. Species richness estimators that 
predict true species richness more accurately are useful 
tools that help to correct this bias. We calculated esti¬ 
mated species richness using the non-parametric ACE 
estimator. Species richness was not estimated for forest 
margin sites due to the small number of spatial replicates 
(n = 2 sites). Species accumulation curves (Mao Tao) 
were calculated to show the number of observed species 
per habitat type dependent on sampling effort. Estima¬ 
tors and species accumulation curves were calculated 
with the Estimates 7.0 package (COLWELL 2006). 

To identify habitat characteristics influencing 
species richness, we measured three habitat parameters 
and ran a multiple ridge-regression model to correct for 
multi-collinearity between explanatory variables. Habi¬ 
tat parameters were canopy closure ‘cc’ (as surrogate for 
shade), slope ‘si’ (as surrogate for flow velocity) and 
stream width ‘sw’ (we decided to measure above bank- 
full, as this factor integrates over a longer period than 
single measurements of water line). Data for slope were 
arcsine transformed before analysis to give a normal dis¬ 
tribution. Several other parameters (dissolved oxygen, 
pH, conductivity) measured during prior studies showed 
little variance among streams (ph ~ 8, dissolved oxygen 
~ 8 mg/L) (Tschelaut 2005, Tschelaut et al. 2005) 
and were therefore excluded from the model. Statistical 
analyses were performed with Statistica 7.0 (StatSoft 
2005) and Primer v5 (Clarke & GORLEY 2001). 

Similarity 

We used Bray-Curtis similarity indices and non- 
parametric multidimensional scaling (NMDS) to visu¬ 
alise relationships of species composition between sam¬ 
pling sites. Sample data were square root transformed 
before analysis. Whether habitat groups had an effect on 
species composition was tested by a one-way analysis of 
similarities (ANOSIM) with Bray-Curtis similarities as 



dependent and the two habitat types f and ca as inde¬ 
pendent factor. Forest margin sites were excluded from 
the ANOSIM because of the small number of spatial 
replicates (n = 2). To examine the potential influence of 
spatial distribution of sites on species composition, we 
performed a Spearman matrix rank correlation for Bray- 
Curtis similarities and distances between sites. 

To test the effects of habitat variables (cc, si, sw) on 
species composition after exclusion of the influence of 
spatial arrangement of sites, we conducted partial Mantel 
tests, using zt software (BONNET & VAN DE PEER 2002). 

Similarity percentage analysis (SIMPER) was used 
to quantify species’ contribution to similarity within a 
habitat and dissimilarity between the three groups of 
habitat types. This method can identify species repre¬ 
sentative for one type of habitat e.g. typical ‘forest 
species’. 

Numbers of individuals and relative species richness 
of widespread and range-restricted species were com¬ 
pared between habitat types using t-tests. Range-re¬ 
stricted (restricted) species were defined as distributed 
only in Central America. Species distributed in Central 
America and additional regions were classified as wide¬ 
spread. Information on the distribution of species was 
extracted from the book ‘Dragonflies and damselflies of 
Middle America and the Caribbean’ (ESQUIVEL 2006). 

Results 

Species richness 

We recorded 902 individuals of adult dragonflies 
representing 32 species belonging to 20 genera and 7 
families (5 Zygoptera, 2 Anisoptera, see Table 1). 

More species were observed in the cultivated land¬ 
scape (18) than at forest sites (11) (t-test: N= 4, t = 4.6, 
df = 3, p = 0.02; Fig. 2) which was also proved by the esti¬ 
mated species richness. Species accumulation curves show 
highest species richness in forest margin habitats, al¬ 
though number of observation units was low (n = 4). Ad¬ 
ditional sampling effort may lead to a higher number of 
species in fm habitats than in other habitat types (Fig. 3). 
Estimated species richness indicates that additional 
sampling effort would add species to the respective 
species inventories (Table 1). However, completeness of 
species inventories was high for f (92%) and ca (90%). 
Estimation of total species richness (53 species) provid¬ 
ed an appraisal of the species inventory that can be ex¬ 
pected in the whole study area representing 19% of 286 
species known for Costa Rica (RAMIREZ et al. 2000). 

The total number of observed individuals was high¬ 
er in cultivated area (638) than in forest (62). Number 
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of observed individuals per site was also higher in culti¬ 
vated area (mean = 159.5, sd = 63.4) than in forest 
(mean = 15.5, sd = 3.9) (t-test: N = 4, t = 4.7, df = 3, 
p < 0.02). A complete list of observations is given in 
Table 2. 

In the multiple ridge-regression model (lambda = 

0.1, multiple r = 0.94, r 2 corr = 0.82, F 3>6 = 17.7, p < 0.004) 

canopy closure, as a surrogate for shade, best accounted 
for variance in species richness (cc: r part = -0.77, 
p = 0.03). Slope and stream width were not significant¬ 
ly related to species richness (sw: r part = 0.62, p = 0.1; 
sb r part = -0.46, p = 0.3, see Fig. 4). 

Similarity 

The NMDS plot based on Bray-Curtis similarities 
showed that forest and cultivated area sites form two 
distinct clusters with forest margin sites situated in be¬ 
tween (Fig. 5). Species composition differed significant¬ 
ly between forest and cultivated area (ANOSIM: global 
R = 0.70, p = 0.003; pairwise tests for cultivated area 
and forest: R = 0.99, p = 0.03). 

Spearman matrix rank correlation with Bray-Curtis 
similarities and geographic distance (Rho = 0.39, 
p = 0.01) outlined the ascending dissimilarity with in¬ 
creasing distance of sites. Partial Mantel tests proved 
that habitat parameters cc and si significantly influ¬ 
enced similarities of dragonfly assemblages (corrected 
for spatial distribution of sites), (cc: r = -0.62, p < 0.001, 
sw: r = -0.02, p = 0.45, and si: r = -0.44, p = 0.001) 
(Fig. 6). 

SIMPER analysis revealed average similarities with¬ 
in groups for forest (39%), cultivated area (36%) and 
forest margin (all similarities were zero). Species con¬ 
tributing most to similarity of habitats were Heteragrion 
erythrogastrum (with 53% contribution to similarity of 
forest), Argia sp. 1 (22%), Hetaerina occisa (15%) for 



Fig. 3: Species 
accumulation curves 
(Mao Tao) for habitat 
types forest margin 
(fm), forest (f) and 
cultivated area (ca). 
Dotted and dashed 
lines without objects 
show 95% confidence 
intervals. 


Table 1: Recorded and estimated (ACE) species for forest (f), cultivated area 
(ca), forest margin sites (fm) and the total survey area. 


Habitat 

Species 

observed 

Estimated 

species 

Coverage% 

Families 

Genera 

f 

11 

12 

92 

5 

7 

ca 

18 

20 

90 

4 

13 

fm 

17 

- 

- 

5 

11 

total 

32 

53 

60 

7 

20 


y = 10.35 -0.07x 
= -0.77, p < 0.03 
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Fig. 4: Habitat 
parameters used in 
the multiple ridge 
regression model to 
account for variance 
in species richness 
between sites, (a) 
Canopy closure was 
found to be a 
relevant predictor for 
species richness, (b) 
Stream width was not 
significantly 
correlated with 
species richness but 
showed a trend to 
more diverse 
assemblages in 
broader streams, (c) 
Correlation between 
slope and species 
richness was not 
significant. 
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Table 2: List of species and number of individuals observed at sampling sites. 


Sites 

Family 

Species 

fl 

f2 

f3 

f4 

fml 

fm2 

cal 

ca2 

ca3 

ca4 

Calopterygidae 

Hetaerina fuscoguttata (Selys 1878) 

- 

- 

- 

1 

- 

7 

6 

20 

20 

40 


Hetaerina occisa (Hagen 1853) 

4 

2 

2 

- 

- 

42 

33 

- 

65 

- 


Hetaerina vulnerata (Hagen 1853) 

- 

- 

- 

- 

- 

- 

10 

16 

40 

5 

Coenagrionidae 

Acanthagrion (Selys 1876) sp. 

- 

- 

- 

- 

- 

- 

- 

- 

- 

1 


Argia cupraurea (Calvert 1902) 

3 

3 

- 

- 

- 

- 

- 

- 

4 

2 


Argia (Rambur 1842) sp. 1 

2 

5 

6 

- 

6 

- 

- 

- 

- 

- 


Argia (Rambur 1842) sp. 2 

- 

- 

- 

2 

- 

4 

- 

- 

- 

- 


Argia (Rambur 1842) sp. 3 

- 

- 

- 

- 

- 

1 

24 

39 

60 

4 


Argia (Rambur 1842) sp. 4 

- 

- 

- 

- 

- 

- 

- 

- 

- 

1 


Enallagma novaehispaniae (Calvert 1907) 

- 

- 

- 

- 

- 

- 

3 

- 

1 

11 


Ischnura capreola (Hagen 1861) 

- 

- 

- 

- 

- 

- 

- 

3 

- 

- 

Megapodagrionidae 

Heteragrion erythrogastrum (Selys 1886) 

6 

4 

8 

2 

- 

4 

- 

- 

- 

- 

Polythoridae 

Cora semiopaca (Selys 1878) 

- 

- 

- 

5 

- 

- 

- 

- 

- 

- 

Protoneuridae 

Neoneura esthera (Williamson 1917) 

- 

- 

- 

- 

- 

- 

- 

- 

- 

1 


Protoneura amatoria (Calvert 1907) 

- 

- 

- 

- 

- 

- 

- 

5 

2 

- 

Gomphidae 

Aphylla obscura (Kirby 1899) 

- 

- 

- 

- 

1 

- 

- 

- 

- 

- 

Libellulidae 

Cannaphila insularis (Kirby 1889) 

- 

2 

- 

- 

- 

- 

- 

- 

- 

- 


Cannaphila vibex (Hagen 1861) 

- 

1 

- 

- 

- 

- 

- 

- 

- 

- 


Dythemis multipunctata (Kirby 1894) 

- 

- 

- 

- 

- 

3 

- 

- 

- 

- 


Dythemis sterilis (Hagen 1861) 

- 

- 

- 

- 

- 

1 

1 

12 

11 

26 


Erythrodiplax connata (Burmeister 1839) 

- 

- 

- 

- 

1 

- 

- 

- 

- 

- 


Erythrodiplax fervida (Erichson 1848) 

- 

2 

- 

- 

1 

- 

- 

- 

- 

- 


Erythrodiplax fusca (Rambur 1842) 

- 

- 

- 

- 

100 

- 

3 

106 

7 

20 j 


Erythrodiplax kimminsi (Borror 1942) 

- 

- 

- 

- 

- 

- 

- 

2 

- 

- 


Libellula foliata (Kirby 1889) 

- 

- 

- 

- 

7 

- 

- 

- 

- 

- 


Macrothemis imitans (Karsch 1890) 

- 

- 

- 

- 

- 

- 

- 

- 

- 

2 


Macrothemis (Hagen 1868) sp. 2 

1 

- 

1 

- 

- 

- 

- 

- 

- 

- 


Nephepeltia phryne (Perty 1834) 

- 

- 

- 

- 

2 

- 

- 

- 

- 

- 


Orthemis ferruginea (Fabricius 1775) 

- 

- 

- 

- 

14 

- 

- 

3 

- 

22 


Perithemis mooma (Kirby 1889) 

- 

- 

- 

- 

- 

1 

- 

5 

- 

- 


Rhodopygia (Kirby 1889) sp. 

- 

- 

- 

- 

7 

- 

- 

- 

- 

- 


Uracis imbuta (Burmeister 1839) 

- 

- 

- 

- 

- 

- 

1 

1 

- 

- 


forest; and Argia sp. 3 (28%), Hetaerina fuscoguttata 
(22%) and Hetaerina vulnerata (15%) for cultivated 
area. 

Average dissimilarity for forest & cultivated area 
was 98% with the most dissimilar species being Argia sp. 
3, Erythrodiplax fusca and Hetaerina occisa. Average dis¬ 
similarity for forest & forest margin was 90% with the 
most dissimilar species being Erythrodiplax fusca , Hetae¬ 
rina occisa and Orthemis ferruginea . Average dissimilarity 
for cultivated area & forest margin was 75% with the 
most dissimilar species being Erythrodiplax fusca , Argia 
sp. 3 and Hetaerina occisa . 

The importance of widespread species increased to¬ 
wards more disturbed habitats (Fig. 7), as mean number of 
individuals and proportion of widespread species were sig¬ 
nificantly higher at sites in cultivated area than in forest 
(t-test; number of species: N = 4, t = 21, df = 3, p < 0.001; 
number of individuals: N = 4, t = 4.6, df = 3, p = 0.02). 


Discussion 

Species richness 

The findings of our study parallel the results of com¬ 
parable research on this issue. Our survey sites located in 
cultivated areas were significantly richer in species and 
had a higher number of individuals overall than forest 
sites. Estimates indicate that a high completeness of 
species inventories was achieved in both habitat types. 
Therefore, we expect a higher sampling effort not to af¬ 
fect the general conclusions drawn from our results. Esti¬ 
mation of total number of species provides an initial as¬ 
sessment of species richness for the area around La Gam- 
ba but supplementary sampling during the rainy season 
will probably improve this relatively low estimate. 

An increase of species richness from primary to¬ 
wards slightly disturbed sites has been found in various 
studies from the subtropics (SAMWAYS & STEYTLER 
1996) and tropics (CLAUSNITZER 2003), and has been 
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explained by the increased habitat diversity created by 
disturbance. Shade and exposed macrophytes are men¬ 
tioned to be the most important environmental param¬ 
eters accounting for variable species composition. Few¬ 
er species are resident in shaded forests than in habitats 
with a light/shade mosaic or open sunny locations. In 
La Gamba we found canopy closure to be the relevant 
influence on species richness. Canopy closure will di¬ 
rectly influence the shading of water surface and banks 
and thus other biotic and abiotic conditions (e.g. tem¬ 
perature and level of dissolved oxygen in the water) as 
important parameters for larvae. Polythorid larvae for 
instance are the only dragonfly larvae in Mesoamerica 
with additional lateral gills on the abdominal segments 
and apparently require well oxygenated water (ESQUIV¬ 
EL 2006). In this family we found Cora semiopaca , a 
species endemic to Costa Rica and Panama, at forest 
site 4, a small, steep stream at the bottom of a shaded 
ravine. 

Disturbed habitats were mostly inhabited by oppor¬ 
tunistic and common species that populate habitats in 
higher densities. Their presence in particular is respon¬ 
sible for the higher species richness at disturbed habi¬ 
tats, while sensitive species vanish. The high species 
number at forest margin habitats may not derive from a 
potential ecotone state of these habitats, combining el¬ 
ements of pristine and altered areas, but is likely to be 
explained by the two different types of surveyed forest 
margin sites. Site fm2 was relatively undisturbed where¬ 
as fml showed a high degree of bank degradation with 
respect to shading trees and was more resembling culti¬ 
vated area sites. 

Although our study has the disadvantages of a low 
sampling effort, a temporal restriction to the dry season 
and an uneven spatial distribution of forest and cultivat¬ 
ed area sites, we do not expect this to affect our general 
conclusions. Most species, however, are known to repro- 



Fig. 5: NMDS plot 
based on Bray-Curtis 
similarities. Sites 
belonging to the 
same habitat category 
(f = forest, 
fm = forest margin, 
ca = cultivated area) 
are connected by 
lines. 
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Fig. 6: Relationship 
between faunal 
similarities and 
canopy cover over 
streams. Similarity of 
species inventory is 
negatively correlated 
to difference in 
canopy closure (cc) 
between sites. 


duce typically during the rainy season showing high 
abundances and some species do not appear as adults 
during dry season. So it is likely that additional species 
not recorded in our study will be found in the survey 
area during the wet season. The potential of other not 
investigated factors influencing the occurrence of 
species will be discussed later on. 



Fig. 7: Proportion of 
species (a) and 
individuals (b) of 
widespread and 
range-restricted 
dragonflies in 
different habitat 
types. Species not 
identified to species 
level are not included 
(3 species for each 
habitat respectively). 
Forest (f): n = 4; forest 
margin (fm): n = 2; 
cultivated area (ca): n 
= 4. Range restricted 
(restricted) species are 
distributed in Central 
America only. 
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Fig. 8: Cultivated area site 4 at Rio Oro, a broad stream with low flow 
velocity. Beyond the adjacent strip of gallery forest the stream is surrounded 
by cattle pastures and paddy fields. 



Fig. 9: Forest site 3 at the stream Quebrada Negra embedded in pristine 
forest. The rapid stream velocity, steep slope and riverine vegetation affect 
inhabiting fauna. 


Similarity 

As expected, we found dissimilar dragonfly commu¬ 
nities at streams in forests, at forest margins and in cul¬ 
tivated areas, representing a gradient in human distur¬ 
bance. Communities in forests and cultivated areas were 
clearly distinct and forest margin sites acted as a link be¬ 
tween the former communities. On the one band, we 
found few species that tended to occur very specifically 
at pristine forest sites (e.g. Heteragrion erythrogastrum ), 
and on the other band many species that were not re¬ 
stricted to distinct habitats and could be encountered 
nearly everywhere (e.g. Hetaerina occisa, Fig. 10). 

Species compositions of three forest sites were very 
similar, most probably because all the sites were situated 
along the same stream. Forest site 4 is characterised by 
different conditions, since there the narrow stream falls 
down a cascade of rocks and then fades in a thin runlet 
into shaded ravine forest. This site is worthy of mention 
for the presence of Cora semiopaca (Fig. 11). SIMPER 
revealed the species responsible for these differences. 
Heteragrion erythrogastrum was frequently found at forest 
streams and forest margin site 2. This species is reported 
to prefer shaded habitats, where it perches for hours in 
the dark (ESQUIVEL 2006). Although all cultivated area 
sites provided shaded patches and cultivated area 1 had 
relatively dense canopy cover, Heteragrion erythrogas¬ 
trum was completely absent from these habitats. It is un¬ 
known to what extent fragmentation and isolation from 
closed forest prevents the colonisation of such streams 
running through small patches of gallery forest by 
species such as Heteragrion erythrogastrum . 

The dragonfly communities of cultivated areas were 
characterised by a significantly higher proportion than 
in forest communities of species with a wide distribution 
(i.e. exceeding the boundaries of Central America). 
Forest sites were inhabited by species specialised to the 
prevailing conditions, but not able to establish popula¬ 
tions in cultivated area, dominated by common and 
widespread species. 

The presence of Erythrodiplax fusca strongly discrim¬ 
inated forest margin habitats from others. We found this 
species in large quantities at forest margin 1. This site, 
although it is situated next to the forest, in some way 
represented the highest degree of degradation of all 
sites, because it was strongly artificially shaped and no 
trees covered its banks. Erythrodiplax fusca is a common 
element of open landscapes on both American conti¬ 
nents (Esquivel 2006). 

Differences between habitats remained on higher 
taxonomic levels since we observed certain families in 
respective habitats exclusively. Polythoridae only oc¬ 
curred at forest streams, Protoneuridae in cultivated area 
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and Gomphidae at forest margin streams. Members of 
Polythoridae and Gomphidae are inhabitants of streams 
and rivers in forests. Especially species of Polythoridae 
are not very common and only species of the genus 
Aphylla (Gomphidae) are reported to live in open areas 
(ESQUIVEL 2006). The differences between dragonfly en¬ 
sembles on species and genus level and the absence of 
some species in the disturbed secondary forest, accompa¬ 
nying forest margin 2 stream, demonstrate the close rela¬ 
tion of some species to specific forest conditions. 

The composition of species assemblages and similar¬ 
ities between habitat types was significantly influenced 
by differences in habitat characteristics, clearly overrid¬ 
ing the effects deriving from the spatial arrangement of 
sites. 

Even though we obtained a clear pattern, some con¬ 
cerns have to be mentioned, since some factors not in¬ 
vestigated in this study may have an additional impact 
on dragonfly communities: 

(1) We conducted our survey for a relatively short 
period of four weeks during the dry season only. The 
likelihood of observing further species during the rainy 
season was mentioned above. 

(2) PlCHLER (2006) reported 23 fish species for the 
Quebrada Negra, some of them insectivorous. Fish pre¬ 
dation on dragonfly larvae is one additional factor shap¬ 
ing the structure of larval composition and distribution 
(JOHANSSON et al. 2006). We have no knowledge of the 
fish fauna in cultivated area streams, but the occurrence 
of insectivorous fish is very likely. A different composi¬ 
tion of predatory fish species in cultivated area sites 
may, therefore, have a different effect on the communi¬ 
ty of aquatic insects. 

(3) Forest streams were low order streams (1-2) and 
the Quebrada Negra was a tributary of the Rio La Gam- 
ba where the site cal was situated. Cultivated area 
streams were of higher stream order and therefore differ¬ 
ences between sites maybe not only be related to habi¬ 
tats but the type of investigated stream and river. Culti¬ 
vated area streams were broad and showed low flow ve¬ 
locities as they ran through a flat landscape, whereas 
forest streams were situated in a mountainous land¬ 
scape. Odonate assemblages usually vary along a water¬ 
course (Hawking & New 1999, Chovanec et al. 2004, 
HOFMANN & MASON 2005) in response to changing 
physical and biotic parameters. Stream width and slope 
gave us no significant response when investigating 
species richness variance between habitat types; howev¬ 
er, these parameters at least show clear trends when re¬ 
lated to species richness, as mentioned before. Further¬ 
more, they are linked to the characteristics of current 
water bodies in plains or mountains, and a part of the 



Fig. 10: Hetaerina occisa, male. This species was encountered in all habitat 
types but was most abundant at streams in cultivated area. 



Fig. 11: Cora semiopaca, male. This endemic species typically lives at shaded 
streams in pristine forests. 
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different community members may be elucidated not by 
deforestation but by a lowland/upland dualism. 

This study focused on dragonfly communities of 
streams, thereby representing only a certain proportion 
of the total species ensemble. However, some species 
living away from flowing waters are not encountered in 
open habitats, including dragonflies specialised to using 
phytotelmata exclusively (i.e. water-filled tree holes and 
tanks of epiphytes) for oviposition and larval develop¬ 
ment. Within the study area, we observed individuals of 
M ecistogaster modestus in the forest. This species is a 
member of the South and Central American distributed 
family of Pseudostigmatidae, whose larvae live in tank 
bromeliads (YANOVIAK 2001) and are described as living 
exclusively in forests. Thus, the occurrence of M ecisto¬ 
gaster modestus stresses the importance of forests as habi¬ 
tats for sensitive species. 

Agricultural land use combined with deforestation 
changes important biotic and abiotic parameters of 
riverine habitats, which strongly influence site selection 
of adult dragonflies and the persistence of larvae (shade, 
macrophyte and riparian vegetation) as well as the 
structure of the surrounding landscape in which adults 
may disperse and forage. These factors determine the 
presence of species and shape community composition. 
The pattern of species distribution we found strengthens 
the assumption that streams in agricultural areas do not 
serve as suitable habitats for many sensitive species that 
are typically found in forests, and that deforestation pos¬ 
es a severe threat to these species. When habitats in de¬ 
graded and open landscapes such as agricultural areas 
are unsuitable as reproductive habitats for sensitive for¬ 
est species, it becomes evident that deforestation en¬ 
dangers not only dragonflies but all aquatic biodiversity, 
implying that the conservation of native habitats de¬ 
serves immediate priority and attention. 
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